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ABSTRACT 


The  IBM  Company's  Continuous  System  Modeling  Program 

was  used  to  simulate  the  lateral  and  directional  flight 

control  systems  of  the  F/A-18  aircraft.   The  model  is 

designed  for  use  in  studies  of  high  angle-of -attack  maneuver 

ing  flight  and  is  restricted  to  the  Auto  Flaps  Up  mode  of 

operation.   The  model  accepts  simulated  pilot  stick  and 

rudder  inputs,  air  data  information,  and  rate  gyro,  angle- 

of-attack,  and  acceleration  feedback  signals.   Outputs  are 

differential  stabilator,  differential  leading-edge  and 

trailing  edge  flap,  aileron,  and  rudder  deflections. 

Typical  input  values  are  used  to  validate  the  model, 
*  • 

generating  output  control  surface  deflections  which  corres- 
pond to  those  expected  for  the  F/A-18  aircraft. 
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I.   INTRODUCTION 

In  this  thesis,  work  completed  at  the  Department  of 
Aeronautics  of  the  Naval  Postgraduate  School  as  part  of  that 
Department's  research  into  the  development  of  advanced 
control  concepts  using  digital  electronic  flight  control  in 
U.  S.  Navy  tactical  combat  aircraft  is  described.   In  con- 
current work,  Carter  [Ref.  1],   discusses  in  part  the 
motivation  for  and  the  scope  of  this  research  program. 
Briefly  reviewing  Reference  1,  it  is  desired  to  expand  the 
application  of  digital  flight  control  technology  to  the 
following  areas: 

1.  Investigation  of  active  control  prevention  of  departure 
from  controlled  flight. 

2.  Comparative  testing  of  new  control  law  algorithms. 

3.  Evaluation  of  modern  control  techniques  such  as  optimal 
control,  observers,  and  model  following  control. 

4.  Simulation  of  degraded  flight  conditions  for  combat 
survivability  studies. 

As  discussed  in  Reference  1,  computer  simulation  of  a 
digital  fly-by-wire  aircraft  was  deemed  the  best  initial 
approach  in  studying  the_ topics  listed  above.   The  Navy/ 
McDonnell  Douglas  F/A-18  is  such  an  aircraft,  and  in  fact  is 
the  first  U.  S.  aircraft  in  production  to  use  a  digital 
fly-by-wire  control  system.   In  Reference  1,  the  computer 


modeling  of  the  longitudinal  axis  control  system  of  the 
F/A-18  is  described.   In  this  thesis,  the  model  is  extended 
to  the  lateral- directional  axes  to  complete  the  control 
system  model.   Future  incorporation  of  an  F/A-18  aerodynamics 
simulation  by  Raithel  [Ref.  2]  and  non-linear  equations  of 
motion  will  result  in  a  complete  "flying"  computer  model  of 
the  F/A-18. 

The  computer  model  was  developed  on  the  Naval  Postgrad- 
uate School's  IBM  Model  3033  general  purpose  mainframe 
computer.   The  code  is  written  for  the  IBM  Company's 
Continuous  System  Modeling  Program  (CSMP)  [Ref.  3],  which  is 
a  Fortran  application  program  designed  to  simulate  dynamic 
systems.   The  text,  "A  Guide  To  Using  CSMP"  by  Speckhart  and 
Green  [Ref.  4],  provides  most  of  the  documentation  necessary 
to  use  CSMP.   CSMP  has  several  advantages  over  a  conventional 
Fortran  program  for  this  simulation.   First,  CSMP  includes 
34  built-in  functions  which  serve  to  model  most  of  the  block 
diagrams  encountered  in  control  systems  engineering.   These 
functions  are  analogous  to  Fortran  functions.   Second,  CSMP 
can  accommodate  user-defined  functions,  called  macros,  and 
can  also  call  standard  Fortran  subroutines  and  functions. 
Third,  output  formatting,  either  printed  or  graphical,  is 
handled  by  CSMP  using  just  a  few  statements.   Finally,  the 
time  base  and  numerical  integration  routines  are  provided 
by  CSMP.   One  might  consider  the  primary  disadvantage  of 
CSMP  to  be  lack  of  user  control  over  the  built-in  functions 
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and  numerical  methods.   For  this  simulation  CSMP  proved 
to  be  a  very  effective  tool. 

In  Chapter  II,  an  overview  of  the  F/A-18  flight  control 
system  is  provided  and  the  assumptions  made  and  limitations 
of  the  modeling  process  are  discussed.   In  Chapter  III,  the 
methodology  and  nomenclature  of  the  computer  program  are 
outlined.   In  Chapter  IV,  the  tests  and  results  used  to 
validate  the  model  are  presented.   Conclusions  and  recommen- 
dations are  in  Chapter  V.   Appendix  A  contains  block  diagrams 
of  the  F/A-18  flight  control  system  as  modeled,  including  the 
previous  longitudinal  model  developed  by  Carter  [1].   Changes 
in  nomenclature  and  arrangement  to  Carter's  model  have  been 
made  in  incorporating  it  into  the  current  model;  these  are 
essentially  cosmetic.   Appendix  B  is  the  computer  program 
listing  in  CSMP  of  the  F/A-18  flight  control  system  as 
modeled,  again,  readers  of  Carter's  program  will  note  some 
changes  in  its  form  here  to  suit  the  overall  program. 


II.   F/A-f18  FLIGHT  CONTROL  SYSTEM 

The  F/A-18  flight  control  system  is  described  in  detail 
in  the  McDonnell  Aircraft  Company's  F/A-18A  Flight  Control 
System  Design  Report  [Ref.  5].   In  this  chapter,  basics  of 
the  system  are  discussed  and  portions  used  in  the  model  are 
indicated. 

The  flight  control  system  is  a  digital  fly-by-wire 
electronic  control  augmentation  system  which  uses  a  four- 
channel  parallel  network  of  computers,  electronic  circuit 
elements,  and  associated  wiring.   The  system  is  entirely 
electronic  from  pilot  controls/feedback  sensors  to  the 
control  surface  actuators.   The  actuators  are  redundant 
electrohydraulic  servo  mechanisms  with  the  exception  of  the 
leading-edge  flap  system  which  is  a  hydraulically-powered 
rotary  mechanical  system.   Backup  mechanical  control  of  the 
stabilators  is  available,  and  the  stabilator,  aileron  and 
rudder  surfaces  have  a  backup  analog  Direct  Electric  Link 
in  the  event  of  a  total  digital  computer  failure.   Control 
surfaces  are  stabilators,  ailerons,  dual  rudders,  leading- 
edge  and  trailing-edge  flaps.   The  stabilators,  leading-edge 
and  trailing-edge  flaps  are  capable  of  differential  movement 

Pilot  inputs  are  through  a  conventional  control  stick 
and  rudder  pedal  arrangement.   Closed-loop  stability 
augmentation  is  provided  by  feedback  of  pitch,  roll,  and 
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yaw  rates,  normal  and  lateral  acceleration,  and  angle-of- 
attack  data.   These  feedback  signals  are  gain  scheduled  by 
air  data  and  angle-of- attack  information  to  tailor  the 
feedback  signals  to  the  current  flight  condition.   Cross- 
axis  control  signal  interconnects  are  provided  which 
improve  control  and  feedback  coordination  and  reduce  inertia 
coupling.   Reference  5  should  be  consulted  for  a  more 
in-depth  discussion  of  the  feedback  and  gain  schedule 
design  theory. 

Initially,  the  computer  model  is  to  be  used  for  studies 
involving  only  the  "up  and  away"  flight  mode,  meaning  that 
the  simulated  aircraft  is  in  a  trimmed,  stable  condition  in 
normal  flight  prior  to  the  initiation  of  a  maneuver.   This 
allows  the  model  to  be  reduced'  from  the  full  control  system. 
The  assumptions  which  were  made  in  reducing  the  model  are 
discussed  in  Reference  1,  however,  several  will  be  repeated 
here  : 

1.  The  aircraft  is  in  the  Auto  Flaps  Up  configuration.   The 
control  law  gain  schedules  used  are  designed  for  cruise  and 
combat  maneuvering  flight.   The  leading  and  trailing-edge 
flaps  are  automatically  positioned. 

2.  Inner  loop  control  is  being  used.   This  means  that  the 
pilot  is  the  source  of  control  inputs  through  the  stick  and 
rudders.   Outer  loop  (autopilot)  control  is  not  modeled. 

3.  The  Control  Augmentation  System  is  in  use.   This  is  the 
optimum  situation  indicating  that  the  digital  computers, 
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feedback  sensors,  and  air  data,  acceleration,  and  angle-of- 
attack  measurements  are  all  operating  normally. 
4.   Trim,  external  stores,  speedbrakes,  and  anti-spin 
features  are  not  modeled.   Active  Oscillation  Control  is 
also  left  out,  as  it  is  a  function  of  external  store 
loading . 

The  model  thus  used  is  acceptable  for  simulating  the 
full  range  of  combat  maneuvering  flight  of  which  the  F/A-18 
is  capable.   Take-off  and  approach/landing  phases  are  not 
simulated. 

In  the  McDonnell  Aircraft  Company's  report,  F/A-18 
Flight  Control  Electronic  Set  Control  Laws  [Ref.  6],  the 
software  design  and  nomenclature  of  the  F/A-18  digital 
control  system  is  described.   Those  portions  which  are 
included  in  the  model  are  briefly  reviewed  below,  including 
for  completeness  the  longitudinal  portion  which  was 
discussed  in  Reference  1.   Readers  desiring  in-depth  infor- 
mation of  the  F/A-18  flight  dynamics  and  control  system 
theory  of  operation  should  refer  to  References  5  and  6. 
LONGITUDINAL  AXIS 
The  control  system  consists  of  the  following  paths: 

o   Stick  pitch  input 

o   Pitch  rate  feedback 

o   Normal  accelerometer  feedback 

o   Angle-of -attack  feedback 
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o   Inertial  decoupling  feedback  (roll  rate  *  yaw  rate) 

o   Forward  loop  integrator 
These  paths  are  summed  to  form  the  stabilator  command. 
Except  for  inertial  decoupling  feedback,  all  of  the  paths 
are  used  as  inputs  to  the  forward  loop  integrator  path, 
which  reduces  the  steady-state  difference  between  maneuver 
command  and  aircraft  response  to  zero.   The  signal  paths 
are  gain  scheduled  by  air  data,  angle-of- attack ,  and 
acceleration  functions  which  are  discussed  in  Chapter  III. 
A  digital  notch  filter  is  located  after  the  feedback  summing 
junction  to  attenuate  structural  vibrations.   In  the  Auto 
Flap  Up  mode*  gain- scheduled  angle-of -  attack  information  is 
used  to  provide  maneuvering  flap  commands .- 
LATERAL  AXIS 
The  control  system  consists  of  the  following  paths: 

o   Stick  roll  input 

o   Roll  rate  feedback 

o   Rudder  pedal  interconnect 
The  lateral  command  is  formed  by  the  sum  of  the  signal  paths, 
which  are  gain- scheduled  similarly  to  the  longitudinal  axis 
paths.   The  stick  and  roll  rate  paths  incorporate  digital 
notch  filters  to  attenuate  structural  vibration  inputs. 
Differential  stabilator,  and  leading-edge  and  trailing  edge 
flap  commands  are  separately  gain-scheduled  and  directed  to 
their  respective  paths  in  the  longitudinal  axis  control 
system. 
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DIRECTIONAL  AXIS 

The  control  system  consists  of  the  following  paths: 

o   Rudder  pedal  yaw  input 

o  Yaw  rate  feedback 

o   Lateral  acceleration  feedback 

o   Rolling  surface  interconnect 

o   Inertial  decoupling  feedback  (pitch  rate  *  roll  rate) 
The  rudder  command  is  the  sum  of  structural-filtered,  gain 
scheduled  feedback  signals,  and  gain-scheduled  rudder  pedal 
and  rolling  surface  interconnect  signals. 

Actuator  data  from  Reference  5  was  not  expressed  in 
terms  of  damping  and  natural  frequency,  therefore,  data 
that  was  available  on  an  actuator  very  much  like  the  F/A-18 
stabilator  actuator  was  used  as  a  guide.   It  was  decided  to 
use  a  second  order  model  for  all  of  the  control  surface 
actuators,  and  to  use  the  same  damping  ratio  and  natural 
frequency,  0.7  and  40Hz,  respectively,  until  the  data 
specified  in  that  manner  was  obtained  for  the  other 
actuators . 

Digital  filters  other  than  the  structural  filters 
already  mentioned  include  lag,  lead-lag,  and  integrators. 
The  constants  for  all  of  the  filters  are  listed  in 
Chapter  16  of  Reference  5.   Aliasing  filters  were  treated 
as  analog  filters  (s-domain)  as  they  always  occur  prior  to 
A/D  converters.   Digital  filters  used  the  z-domain 
representation. 
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In  Chapter  III,  the  modeling  methodology  by  which  the 
F/A-18  digital  control  system  as  described  briefly  above 
was  coded  for  computer  study  using  CSMP  is  discussed. 
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Ill .   PROGRAM  METHODOLOGY 

The  primary  source  of  information  used  in  making  the 
computer  model  was  the  McDonnell  Aircraft  Company  F/A-18 
Flight  Control  System  Design  Report  [Ref.  5],   Figures 
16.1,  16.2,  and  16.3  of  Reference  5  are  the  block  diagrams 
of  the  longitudinal,  lateral,  and  directional  control 
systems,  respectively.   Chapter  16  of  Reference  5  contains 
descriptions  of  system  operation,  control  law  algorithms, 
and  digital  filter  specifications.   Since  there  have  been 
several  versions  of  the  control  laws  to  date,  it  should 
be  noted  that  the  version  used  in  this  program  is  OPV  8.2.1, 
current  as  of  31  August  1982. 

Program  modeling  began  by  reducing  the  control  system 
block  diagrams  (Figures  16.1,  16.2,  16.3  of  Reference  5) 
to  the  forms  desired  for  research.   This  was  done  by 
applying  the  assumptions  listed  in  Chapter  II  and  in 
Reference  1.   Next,  scheme  of  labeling  the  control  system 
paths  was  developed.   Readers  of  Reference  1  should  note 
that  the  labeling  system  used  there  has  been  changed  in 
this  program  to  accommodate  the  more  general  nature  of  the 
full  three-axis  model.   The  paths  which  are  modeled  are: 

1.  Pilot  inputs:  pitch,  roll,  and  yaw 

2.  Angle-of -attack  feedback 

3.  Rate  gyro  feedbacks:  pitch,  roll,  and  yaw 
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4.  Normal  accelerometer  feedback 

5.  Lateral  accelerometer  feedback 

6.  Stabilators ,  symmetrical  and  differential 

7.  Ailerons 

8.  Rudders 

9.  Leading-edge  flaps,  symmetrical  and  differential 
10.   Trailing-edge  flaps,  symmetrical  and  differential 

The  labels  and  nomenclature  for  each  path  are  listed  below 

PILOT  INPUTS 

PP       pilot  pitch  path 

PR       pilot  roll  path 

PY       pilot  yaw  path 

ANGLE -OF -ATTACK  FEEDBACK 

AA       angle-of- attack  path 

ALPHAT    computed  angle-of -attack 

RATE  GYRO  FEEDBACKS 

PG       pitch  rate  gyro  path 

RG       roll  rate  gyro  path 

YG       yaw  rate  gyro  path 

CSAOAT    approx.  cos (ALPHAT) 

SNA OAT    approx.  sin (ALPHAT) 

P2A,B,C   filter  P2  arguments 

Y3A,B,C   filter  Y3  arguments 

NORMAL  ACCELEROMETER  FEEDBACK 

NZ       normal  accelerometer  path 

NZA      nz  for  gain  schedule  use 
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NZAF  filtered  incremental  nz 

NZAR1  nza-based  roll  rate  gyro  parameter 

P5A,B,C  filter  P5  arguments 

LATERAL  ACCELEROMETER  FEEDBACK 

NY  lateral  accelerometer  path 

STABILATOR  PATH 

ST  main  stabilator  path 

DS  differential  stabilator  path 

ST2R  main  to  differential  stabilator  path  signal 

RST  right  stabilator  path 

RSTDEF  right  stabilator  deflection 

LST  left  stabilator  path 

LSTDEF  left  stabilator  deflection 

AILERON  PATH 

AL  main  aileron  path 

YV1R  rudder  to  roll  crossfeed  path 

RAL  right  aileron  path 

RALDEF  right  aileron  deflection 

LAL  left  aileron  path 

LALDEF  left  aileron  deflection 

RUDDER  PATH 

RD  main  rudder  path 

RSR  roll  surface  to  rudder  interconnect  path 

LRD  left  rudder  path 

LRDDEF  left  rudder  deflection 

Y5A,B,C  filter  Y5  arguments 
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LEADING-EDGE  FLAP  PATH 

LE       main  leading-edge  flan  path 

CO  i.     r 

DLE      differential  leading-edge  flap  path 

RLE      right  leading-edge  flap  path 

RLEDEF    right  leading-edge  flap  deflection 

LLE      left  leading-edge  flap  path 

LLEDEF    left  leading-edge  flap  deflection 

TRAILING-EDGE  FLAP  PATH 

TE       main  trailing- edge  flap  path 

DTE      differential  trailing  edge-flap  path 

RTE      right  trailing-edge  flap  path 

RTEDEF    right  trailing-edge  flap  deflection 

LTE      left  trailing-edge  flap  path 

LTEDEF    left  trailing-edge  flap  deflection 

Nomenclature  from  the  McDonnell  Aircraft  Company  literature 

[Refs.  5  and  6]  has  been  included  where  it  was  deemed 

helpful  for  reference.   This  primarily  includes  such  terms 

as  PK_,  RK_,  YK_,  PV_,  RV_,  and  YV_;  these  are  constants  or 

variables  which  occur  at  significant  points.   The  terms 

QC  (dynamic  pressure)  and  PS  (static  pressure)  are  also 

from  the  McDonnell  Douglas  Aircraft  Company  literature. 

Carter  [Ref.  1]  lists  and  defines  some  additional  terms 

which  have  been  used  in  this  program,  including  filter  and 

constant  nomenclature.   The  remainder  of  the  terminology 

was  developed  specifically  for  this  program.   Block  diagrams 

incorporating  the  simplying  assumptions  mentioned 
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previously  and  using  the  above  notation  are  included  in 
Appendix  A.   These  particular  block  diagrams  are  intended 
only  to  aid  in  identifying  the  various  terms  with  their 
location  in  the  diagrams  and  leave  out  much  of  the  detail 
found  in  Figures   16.1,  .2,  .3  of  Reference  5. 

Carter  [Ref.  1,  Chapter  III]  gives  an  excellent  des- 
cription of  the  methodology  of  using  CSMP  to  model  the 
block  diagram  elements  of  Figures  16.1,  16.2,  and  16.3 
of  Reference  5.   Some  changes  have  been  made  to  those 
procedures,  however: 

1.   The  control  law  gain  schedules,  referred  to 

singularly  as  Function  F ,  have  been  removed  from  nosort 

sections.   Instead,  all'of  the  gain  schedules  have  been 

coded  as  Fortran  functions  and  called  as  required.   This 

step  has  removed  over  300  lines  of  code  from  the  CSMP 

program,  since  Fortran  subroutines  or  functions  do  not 

count  against  the  allowable  number  of  CSMP  statements. 

Twenty-four  functions  have  been  added  for  the  lateral  and 

directional  axes ,  in  addition  to  13  for  the  longitudinal 

axis  discussed  in  Reference  1: 

FUNCTION  DESCRIPTION 

F4       Roll  Rate  Feedback  Gain  Schedule  (QC,PS) 

F6       Differential  Stabilator  Gain  Schedule  (RI,PS,RV11) 

F7       Lateral  Command  Gain  Schedule  (QC,PS) 

F10       Rudder  Command  Gain  (QC) 

F13      Lateral  Command  Gain  Schedule   ( RI , PS , STORES) 
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F17      Rudder  Pedal  Command  Gain  Increment  (AOA) 

F30       RSRI  Gain  Schedule  (QC,PS) 

F31       Differential  T.E.  Flap  Gain  Schedule  (PJ,PS) 

F34       Differential  T.E.  Flap  Gain  Schedule  (AOA) 

F35      Lateral  Forward  Loop  Gain  Schedule  (AOA) 

F36      Aileron  Gain  Schedule  (QC,PS,RI) 

F38      RSRI  Gain  Schedule  (AOA,RI,PS) 

F39      Rudder-Roll  Interconnect  Gain  Schedule  (AOA) 

F41      Roll  Surface  Limit  Schedule  (A0A,RI) 

F42       RSRI  Nonlinear  Gradient 

F45       Directional  Forward  Loop  Gain  Schedule  (QC,PS) 

F90       Lateral  Acceleration  Feedback  Gain  Schedule  (RI,PS) 

F93      Differential  L.E.  Flap  Gain  Schedule  (RI,PS,NZAF) 

F96      Yaw  Rate  Gain  Schedule  (QC,PS) 

F101      Differential  Stabilator  Load  Alleviation  (RI,PS, 

NZAF) 
F108      Directional  Inertial  Gain  Schedule  (QC) 
F112      Lateral  Acceleration  Gain  (RI) 
F113      Lateral  Acceleration  Gain  (AOA) 
F114      Rudder  Pedal  Command  Gain  Increment  (RI) 
Algorithms  for  the  implementaion  of  the  functions  are  found 
on  Chapter  16  of  Reference  5. 

2.   Frequency  averagers  and  rate  limiters  are  modeled 
using  CSMP  macro  statements.   Thus,  all  nosort  sections  have 
been  eliminated  from  the  current  3-axis  program. 
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3.  Cross-axis  signal  paths  for  differential  stabilators, 
leading-edge  and  trailing;- edge  flaps  have  been  added. 
Rolling-surface- to- rudder  interconnect  (RSRI),  rudder-to- 
roll  crossfeed,  and  inter-axis  distribution  of  rate  gyro 
feedback  signals  have  been  incorporated. 

4.  As  noted  in  Reference  1,  in  the  longitudinal -only 
simulation,  the  stabilators  and  flaps  were  modeled  for  one 
direction  of  motion  only;  this  constraint  has  been  removed 
in  the  current  3-axis  program. 

Descriptions  of  the  algorithms  used  to  implement  A/D 
and  D/A  converters,  digital  filters,  aliasing  filters, 
limit  functions ,  and  actuator  servomechanisms  are  contained 
in  Reference  1.   They  are  used  in  the  current  program 
without  modification. 
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IV.   MODEL  VALIDATION 

The  process  of  validating  the  computer  model  of  the 
F/A-18  flight  control  system  was  carried  out  in  two  major 
steps  : 

1.  Individual  block  diagram  elements  were  tested  using 
probable  ranges  of  input  values.   Digital  filter  macros, 
rate  limiters,  frequency  averagers,  and  gain  schedule 
functions  are  included  here.   Carter  [Ref.  1]  describes 

the  testing  of  all  of  the  block  diagram  elements  mentioned 
with  the  exception  of  the  gain  schedule  functions  relevant 
to  the  lateral  and  directional  axes,  which  were  individually 
verified  by  this  researcher  exactly  as  were  the  gain 
schedule  functions  used  in  Reference  1. 

2.  The  entire  model  was  assembled  and  subjected  to 
inputs  from  the  stick,  rudder,  and  angle-of- attack , 
acceleration,  and  feedback  sensors.   The  goal  of  this  step 
was  to  verify  correct  direction  of  motion  of  control 
surfaces  in  response  to  unambigous  inputs. 

Inasmuch  as  Reference  1  has  already  discussed  the  procedures 
and  results  relevant  to  Step  1,  with  the  exceptions  noted, 
the  current  discussion  will  involve  Step  2.   Further, 
since  Reference  1  contains  the  validation  results  for 
stabilator,  leading-edge  and  trailing-edge  flaps  due  to 
pitch  inputs,  here  the  control  surface  motion  due  to  roll 
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and  yaw  inputs  will  be  of  primary  concern.   The  source 
of  information  used  as  a  reference  in  comparing  model 
performance  to  design  specification  was  the  McDonnell 
Aircraft  Company  F/A-18  Flight  Control  System  Design 
Report  [Ref.  5]. 

The  computer  model  requires  the  following  inputs : 

1.  Static  pressure  (PS)  in  lb/sq.ft. 

2.  Dynamic  pressure  (QC)  in  lb/sq.ft. 

3.  Pilot  pitch  input  (PP1)  in  lb.  +  =  aft  stick. 

4.  Pilot  roll  input  (PR1)  in  lb,  +  =  right  stick. 

5.  Pilot  yaw  input  (PY1)  in  lb ,  +  =  right  rudder. 

6.  Pitch  rate  gyro  feedback  (PG1)  in  deg/sec, 
+  =  nose-up. 

7.  Roll  rate  gyro  feedback  (RG1)  in  deg/sec, 
+  =  right  roll. 

8.  Yaw  rate  gyro  feedback  (YG1)  in  deg/sec, 
+  =  right  yaw. 

9.  Angle-of -attack  (AA1)  in  deg. 

10.  Normal  acceleration  (NZ1)  in  "g",  +  =  nose-up  motion 

11.  Lateral  acceleration  (NY1)  in  "g",  +  =  nose-right 
motion. 

Depending  on  the  need,  inputs  can  be  programmed  in  CSMP  as 
step  or  ramp  functions  (see  Reference  4) ,  or  using  any 
standard  Fortran  function  such  as  SIN  or  EXP,  for  example. 
Combinations  thereof  are  acceptable,  as  well.   An  example 
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of  a  test  roll  input  consisting  of  a  61b  right  stick  input 
for  two  seconds  followed  by  a  61b  left  stick  input  as 
written  in  CSMP  would  be: 

PR1  =  6.0*STEP(0.0)  -  12.0*STEP(2  .0) 

For  this  purposes  of  the  tests  described  in  this  thesis, 
PS  and  QC  are  held  constant  through  test  maneuvers,  though 
there  exists  no  such  constraint  in  the  model.   The  effects 
of  varying  PS  and  QC  were  verified  to  be  correct  when  each 
gain  schedule  function  was  compared  to  the  graphical  data 
in  Chapter  16  of  Reference  5.   Rate  gyro  feedback,  angle- 
of-attack,  and  acceleration  inputs  may  currently  be  pro- 
grammed as  described  for  pilot  inputs.   However,  accurate 
representation  of  these  latter  three  types  of  inputs  will 
not  be  possible  until  the  F/A-18  aerodynamic  build-up  by 
Raithel  [Ref.  2]  and  the  non-linear  equations  of  motion 
are  incorporated  into  the  program.   As  that  time,  those 
inputs  will  be  determined  by  the  program  and  will  not  be 
explicity  stated.   For  the  tests  described  here,  rate  gyro 
feedback,  angle- of -attack ,  and  acceleration  data,  when 
necessary  to  simulate  aircraft  responses,  will  be  provided 
by  "best  guess"  estimation  of  those  data,  given  the  pilot 
input . 

Figures  presented  in  this  chapter  are  based  on  a  time 
scale  of  four  seconds  per  maneuver.   Control  surface 
deflections  are  in  degrees.   Static  pressure  and  dynamic 
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pressure  inputs  are  listed  on  each  figure  in  terms  of 
altitude  and  Mach  number.   Standard  day  conditions  are 
assumed  in  all  cases.   Angle-of -attack  is  occasionally 
listed  where  it  plays  a  significant  role  in  shaping  the 
control  system  response. 

Figure  4.1.  depicts  differential  stabilator  response  to 
a  +/-  6.01b  lateral  stick  input.   The  shift  from  right  to 
left  stick  takes  place  at  the  2.0  second  mark.   The  per- 
tinent gain  schedule  functions  are  Functions  6  and  101. 
For  the  given  flight  conditions  and  pilot  input,  the 
stabilator  response  is  very  close  to  maximum.   The  response 
is  decreased  at  higher  dynamic  pressures. 
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Figure  4.1 


Differential   Stabilator   Response    to    Roll 
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Figure  4.2.  shows  differential  stabilator  response  to 
the  same  input  as  for  Figure  4.1.,  but  with  roll  rate  gyro 
feedback  added  as  a  ramp  type  input.   This  causes  the  damp- 
ing effect  on  the  stabilator  motion.   The  overshoot  caused 
when  the  input  signals  are  reversed  in  sign  is  to  some 
extent  believed  to  be  an  artifact  of  the  programming  of 
the  roll  rate  feedback  input,  which  only  approximates  an 
actual  feedback  input.   This  phenomena  can  be  noticed  on 
some  of  the  following  graphs,  as  well. 

Figure  4.3.  shows  differential  stabilator  motion  due  to 
a  rudder  input.   This  demonstrates  the  rudder- to- roll 
surface  CAS  interconnect,  which  is  gain-scheduled  by 
Function  39  using  angle-of -attack .   A  ramp  input  was  used 
to  simulate  the  angle-of -  attack  signal. 

Aileron  motion  due  to  a  +  /-  6.01b  roll  input  is 
depicted  in  Figure  4.4.   Aileron  gain  schedules  are 
Functions  35  and  36,  which  accept  angle-of -  attack  and  air 
data  inputs,  respectively.   Both  functions  decrease  gain 
as  their  input  values  increase:  in  the  case  of  increasing 
angle-of -attack ,  this  is  to  reduce  sideslip;  the  air  data 
gain  schedule  lessens  aileron  response  at  high  dynamic 
pressures  where  aeroelastic  effects  would  otherwise  cause 
aileron  reversal. 
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Figure  4.2.   Differential  Stabilator  Response  with 
Roll  Rate  Feedback 

When  roll  rate  feedback  is  added  to  the  aileron  roll 
motion, it  is  seen  that  a  damping  effect  is  present,  as  it 
should  be.   This  is  shown  in  Figure  4.5.   Roll  rate  feedback 
is  gain- scheduled  by  Function  4,  which  acts  to  decrease 
feedback  gain  as  dynamic  pressure  increases.   The  overall 
lateral  control  system  gain  is  shaped  additionally  by 
Functions  7  and  13,  using  air  data  inputs.   These  gain- 
schedules  do  not  easily  lend  themselves  to  straight  -  forward 
explanations  of  purpose,  indeed,  there  are  certainly  multi- 
ple purposes  for  the  design  of  these  functions.   Reference  5 
is  more  explicit  as  concerns  the  gain-schedule  design 
philosophy . 

28 


"1 

I 

PI 
»1 


UJ 

o 


%DIFF.STR3.     RESPONSE    TC    ""20    LB    TAW 
'    M=0.5Q.H=5QOOFT.flCR=9.6+l. 75    DEG/SEC 


o.oo    a.  no    o.so 


1-20      1.60      2.00      2.H.0      2.BD      3.20      3.60      U.OO 
TIME 


LEGEND 
RSTOE^   2 
LSTOEF   O 


Figure    4.3.       Differential    Stabilator    Response   to   Yaw 
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Figure  4.4.   Aileron  Response  to  Roll 
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Figure  4.5.   Aileron  Response  with  Roll  Rate  Feedback 
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Figure  4.6.   Aileron  Response  to  Yaw 
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The  rudder- to- roll  surface  CAS  interconnect  signal 
causes  an  aileron  deflection  as  angle -of -attack  is 
increased.   This  effect  is  shown  in  Figure  4.6.,  and  is 
analogous  to  that  noted  in  Figure  4.3. 

At  higher  dynamic  pressures,  as  noted  earlier,  aileron 
gain  is  decreased  to  reduce  aeroelastic  effects.   Differen- 
tial leading-edge  and  trailing-edge  flaps  are  incorporated 
to  maintain  an  acceptable  roll  response  under  these  cir- 
cumstances.  Differential  leading-edge  flap  motion  in 
response  to  a  +/-  10.0  lb  roll  input  is  displayed  in 
Figure  4.7.   At  the  indicated  Mach  number,  aileron  response 
is  nil.   Function  93  governs  the  response  as  a  function  of 
air  data  and  normal  acceleration  inputs. 
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Figure  4.7.   Differential  Leading-Edge  Flap  Response  to  Roll 
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Differential  Leading-Edge  Flap  Response  to  Roll 
with  Roll  Rate  Feedback 


The  expected  effect  of  roll  rate  feedback  on  differen- 
tial leading-edge  flap  motion  is  shown  in  Figure  4.8.   The 
damping  is  noticeably  less  than  for  the  differential 
stabilator  and  aileron  damping  observed  previously;  this  is 
due  to  the  inherently  higher  damping  of  aircraft  rolling 
motion  at  the  dynamic  pressures  where  differential  flap 
motion  is  required.   Figures  4.9.  and  4.10.  depict  similar 
responses  for  differential  trailing- edge  flaps,  without  and 
with  roll  rate  feedback,  respectively.   Functions  31  (air 
data)  and  34  (angle-of- attack)  are  the  gain  schedules 
responsible  for  differential  trailing-edge  flap  motion. 
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Figure    4.9.       Differential   Trailing-Edge    Flap    Response    to    Roll 
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Figure    4.10.    Differential   Trailing-Edge    Flap    Response    to    Roll 
with    Roll    Rate    Feedback 
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Figure  4.11.   Rudder  Response  to  Yaw 

Directional  control  system  responses  are  made  by  the 
twin  rudders,  which  do  not  act  differentially  in  the  Auto 
Flap  Up  mode.   Rudder  response  to  a  +/-  20.01b  rudder  pedal 
input  is  shown  in  Figure  4.11.   The  gain  schedules  which 
shape  the  motion  are  Functions  10  (air  data),  71  (air  data), 
and  114  (angle-of -attack) .   Figure  4.12.  shows  rudder  motion 
when  feedback  signals  are  added  to  the  simulation.   The 
damping  effect  is  the  sum  of  stability  axis  yaw  rate  (this 
feedback  signal  is  a  blend  of  yaw  and  roll  rates,  and  angle- 
of-attack  data),  and  lateral  acceleration  feedbacks. 
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Figure  4.12.   Rudder  Response  to  Yaw  with  Feedback 
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Figure  4.13.   Rudder  Response  to  Roll 
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Several  gain  schedules  are  involved  here  in  shaping 
rudder  response  to  the  requirements  of  sideslip  reduction, 
prevention  of  excessive  vertical  tail  loads,  inertia 
coupling  reduction,  and  high  angle-of- attack  maneuvering 
stability.   Reference  5  discusses  these  considerations. 

Rolling  surface-to-rudder  interconnect  signal  operation 
can  be  seen  in  Figure  4.13.   Functions  30  (air  data),  and 
38  (angle-of -attack  and  air  data)  shape  this  response. 

Verification  of  all  of  the  signal  paths  of  the  F/A-18 
flight  control  system  in  the  Auto  Flap  Up  mode  has  thus 
been  accomplished.   The  gain  schedules  (Functions)  were 
each  individually  verified  prior  to  incorporation  into  the 
computer  model . 

The  process  of  model  validation  has  thus  been  to 
observe  model  responses  to  some  relatively  basic  inputs, 
and  then  to  judge  these  responses  as  being  plausible  or  not 
with  respect  to  the  data  provided  in  Reference  5.   The 
computer  model  has  given  logical  results  as  noted  in  this 
chapter,  and  also  for  a  wider  range  of  input  signals  and 
flight  conditions  than  have  been  discussed  here.   The 
conclusion  is  that  the  flight  control  svstem  model  is  valid 
representation  of  the  F/A-18  flight  control  system  within 
the  range  of  the  simplifying  assumptions  and  limitations 
noted  previously  in  Reference  1  and  in  this  thesis. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  computer  simulation  of  the  F/A-18  digital  electronic 
flight  control  system  developed  here  is  suitable  for  further 
use  in  aircraft  control  systems  studies.   With  the  inclusion 
of  F/A-18  aerodynamics  and  equations  of  motion,  it  will  also 
serve  to  examine  simulated  maneuvering  flight  under  condi- 
tions found  near  or  at  the  accepted  flight  envelope.   Com- 
parative studies  may  be  done  of  newer  control  systems 
concepts,  using  the  known  response  as  the  basis  for 
comparison. 

CSMP  has  size  limitations  for  several  internal  para- 
meters, the  most  pertinent  of  which  are  total  number  of 
program  statements  (1900),  and  maximum  number  of  statements 
in  any  given  sort  section  (600)  .   The  current  program  is 
right  at  the  latter  limit  with  599  statements  in  one  sort 
section.   Some  statements  could  be  combined  to  reduce  this 
number;  this  is  not  recommended  until  the  user  is  familiar 
with  the  F/A-18  flight  control  system  and  its  representation 
here.   It  is  recommended  that  the  aerodynamics  and  equations 
of  motion  be  written  as  Fortran  subroutines  and  called  at 
the  beginning  of  the  dynamic  section.   This  would  use  only 
two  lines,  thus  only  two  additional  statements  would  have 
to  he  combined  to  make  the  necessary  room.   There  are  not 
obvious  dividing  lines  where  the  current  program  could  be 
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broken  up  into  more  than  one  sort  section,  however,  users 
experienced  in  CSMP  may  find  this  to  be  possible  when  the 
aerodynamics  and  equations  of  motion  are  incorporated. 
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APPENDIX  A:      MODEL   BLOCK   DIAGRAMS 


PILOT    PITCH    INPUT    PATH 


pilot  otten 
•tic*  forct 


PBl 


■*cn«n1ca1 


DP  2 


•  t'ot  dyna»1c* 


PP3 


aliasing  inter 


Dt>4 


a/O  converter 

(M   nz) 


PP5 


stick  araotant 
(function  20) 


PP« 


Maltar 


PD? 


pv3 


V 
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PILOT  ROLL  INPUT  PATH 


forc«     * 

I 

•/d  conv«rt»r 
(SO  hz) 

I            J 

I 

' 

pri 

pr6 

■wcmnfcai 

Dr««kout 

CMtdCMnO 

\ 

DT2 

pr7 

•tic*  Oyraalc* 

•tick   grtOt»nt 

i 

PT3 

r 

J 

DTt 

IfaftOT 

ftlt«r  R2 

' 

OT4 

. ! 

DT9 

•Ifaatng  ffltcr 

pr5 

0T9 

\x 

40 


PILOT  YAW  INPUT  PATH 


p 

ruOOtr 


»tr   fore*    \ 


Dyl 


trMkout 


Py2 


•  tfasfng  Mltr 


py3 


«/d  eor»v«ct#r 
(40  hx) 


oy« 


aacOMnd 


oy« 


ruooe*"  vat) 
(fincUoo  1) 


0.5-fl7*f 114 


py« 


-B 


oy7 


H»lt«r 


pyt 


py« 

V 
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PITCH    RATE    GYRO   PATH 


pttcn  rttt 
tjyro 


Pfll 


«U««lng   fflttr 


P02 


•/a  cooYirnr 
(•o  nz) 


pg3 


ftu«r  P? 


DO* 


»t«ft«r 


pv4 


pgs 


•a«pl«r   (20  r»z) 


PC* 

V 


pv4« 
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ROLL  RATE  GYRO  PATH 


roll  rtt* 
©to 


rgl 


«1U«fng  ffU»r 


ro2 


' 

' 

•/d  converter 

(40  nz) 

rgS 

' 

' 

ftlttr  «1 

rg4 


ff1t«r   83 


rg5 

I 

flU«r  R4 
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YAW    RATE    GYRO    PATH 


y«*  rat* 
flyro 


ygl 


•  lltvfng  ftitr 


yg2 


•/d  eortverttr 
(40  trr) 


yg3 


yg4 


snao«t 


F~ 


7^ 


<- 


•® 


V7 


•aaoler   (40  hz) 


yg5 


rtte   H«ftar 


ygS 


ftltr  y3 


yd  7 


rtt*  llBlter 


44 


ANGLE-OF-ATTACK    SENSOR   PATH 


■ng)a)  of 
attach 


aal 


altastng  fTUar 


•  •2 


•/d  coov«rt»r 
(40  hz) 


^ 


22.0 


•  aS 


llmfter 


•  •? 


t*17 


pvl 


••■oler   (20  ftz) 


F32A»f37 


pvla 

V 


45 


NORMAL  ACCELEROMETER  SENSOR  PATH 


tcc»!rotti» 


nzl 


•  UMlng  f!U»r 


fU2 


a/d 


(40  nz) 


nz3 


Maltar 


0 


•50.0 


3*0.0 


O 


■& 


nz4 


ft  WW  P5 


ra5 


itattac 


pv? 


IS 


nz7 


••■o'««"  (20  ra) 


rt7» 


V 


6 


mhoIvt   (20  nz) 


nz4« 


filter  P16 


3.5«f»2A 


nrt-f 


V 
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LATERAL  ACCELEROMETER  SENSOR  PATH 


e 


nyl 


•1t««fng  ITUer 


ny2 


•/d  convert! 
(40  hz) 


ny3 


MO 


rg3 

V 


M«pl«r    (40  fi2) 


•s^ 


ny4 


11«1ter 


:X 


ny5 


ftU«r   y7 


L4 


(1. 0/347. 0)**2 


nyt 


Fll?*fll3 


m 


ny7 
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STABILATOR  PATH 


F12«f32A0.5 


l«r  (20  hz) 


t*19 


V 


F12*F32A*F8» 


F12«F32A 


••■pier   (20  nz) 


ss«ol«r   (20  nz) 


t*21 


^ 


^ 


st2 


its 

! 

' 

I 

' 

llatttr 

lf»ft#r 

1 

f 

' 

!    , 

si 

4 

si 

1 

1*22 


•  tl 
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STABILATOR   PATH    (CONT.) 


0.5 


ISI 


-&■ 


«t< 


»ta« 


fUltr  ?9 
(fwd.  trrttgrator) 

•t7  t«prfcit 


*ta 


1! utter 


20   to  40 


V 


•  19 


\y 
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STABILATOR  PATH  (CONT.) 


»t9 

y 

Dvl 

Y     :, 

DV3 


40    to   80 
aver«g«r 


•tl3 


»tu 


Matter 


•t!5 


ftU«r  PI 


•US 


J     «tl7 

Matter 

1 

' 

•t 

IS 
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STABILATOR    PATH    (CONT.) 
DIFFERENTIAL    STABILATOR   SIGNAL    PATH 


11«lt«r   (F4l) 


o*3 


W-«t2r2) 


0*4 


rt\0   H«tt«r 


«t2r2 


0*5 


\y 
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STABILATOR  PATH  (CONT.) 


S118 


1»tl 


ifalttr 


l»t2 


d/«   Conv»rt»r 
(10  hz) 


l«t3 


l#ft   •t«Dt1«tor 
•ctu«tor 


0»5 


rati 


1 luttBT 


rst2 


<3/»  converter 
(•0  trz) 


r*t3 


PtjSTit 

•tt&nator 

actu«tor 


lttorf 


r»tarf 
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AILERON  PATH 


(F4-fl3)*f7 


rv2i 


l»«ft»r 


rvZ 


Dr9 


F31»F34 


rvll 


ft 

(rvll) 


MIN(rkl7*f35.rkl 
7-F101) 


Pkl7 


rv7 


rv7 


1.33 


yvlr 


Unlter 


rvl 


40  to  80 
mv«rig«r 


K 


rvll 
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AILERON   PATH    (CONT.) 


« 

1 

f35*F36 

' 

' 

1  laltM' 

•  12 

X 

■  13 

1 

' 

Itufttr   (F4l) 

_J 

• 

X 

9.5 

• 

4 

< 

' 

o*tOD«nd 

A 

> 

115 


815 


\y 
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AILERON   PATH    (CONT.) 


Tinner 


Tall 


rate  Halter 


Tal2 


a/9  converter 
(»0  rrz) 


T«T3 


Teft  aileron 
actuator 


Tlntter 


rail 


rate  Heller 


ra!2 


oVa  converter 
(80  nz) 


ra13 


ri  cf>t   af  leron 
actuator 


leidef 


raiaef 
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RUDDER   PATH 


^ 


0.5 


F30«f3« 


W 


r«rl 


M»o'r-   (40  r>z) 


rv2 


fTUer  Y5 


r«r3 


•El 


rmrA 


F*2   (rr4) 


0«5 


56 


RUDDER   PATH    (CONT.) 


t>«3 


ny7 


VB* 


saapler   (40  rtz) 


V      V 


»«»oler   (40  nr) 


13 


noe 


yv3d 


Hotter 


IP 


T96 


yv3c 


1!»1t»c 


ftU#c   Yl 


y>3« 


yyJG 


fflt»r   Y2 


yv3 
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RUDDER   PATH    (CONT.) 


ll«(t«r 


lroi 


40  to  80 
■vertper 


lrd2 


d/t  conv*rt«r 
(•0  ftz) 


^r03 


l«ft  rucw«r 
•Ctu«toc 


py« 


F45 

__^7 

r02 

, 

^ 

' 

11a1t«r 

I 

rdi 


1l«ttar 


rrdl 


40  to  80 
•veciger 


rrd? 


d/«   converter 
(80  m) 


•d3 


Mgnt  ruooer 
•Ctuttor 


lrooef 


rrddrf 
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LEADING-EDGE    FLAP   PATH 


•  iDTut 


ftlt»r  PH 


1*2 


«««Ol«r    (20  HZ) 


l»2t 


F27 


1«2tl 


1t«tt«r  (F2») 


l»2t2 


lt«U«r    (P29U) 


1*3 


r«t»   11»U»r 


i»4 


20   to   10 
•v«r«ger 


1*5 


V 
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LEADING-EDGE  FLAP  PATH  ( CONT . ) 
DIFFERENTIAL  LEADING-EDGE  FLAP  PATH 
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LEADING-EDGE    FLAP   PATH    (CONT.) 


n«l 


ltatter 


11«2 


rate   Halter 


U«3 


ltafter 


11* 


<J/«  converter 
(«0  nz) 


lli-i 


left 

lea<Jfnq-#<jg« 

fl«0  actuator 


ie5 

V 


dle3 


Mel 


ltatltr 


rU2 


" 

rate  Halter 

i 

(  rle3 

Ifafter 

r1«4 


d/e  converter 
(80  hz) 


tt^t 


ptgnt 
leedlnc-eoge 
flip  ectuator 


lledef 


rleOef 
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TRAILING-EDGE    FLAP   PATH 


alDTMt 


ftUer  P12 


t»2 


««*oler   (20  nz) 


t»2t 


F24 


t»2tl 


llwtter   (F25) 


t«3 


20   to  SO 
•veragw 
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TRAILING-EDGE  FLAP  PATH  (CONT.) 
DIFFERENTIAL  PATH  INCLUDED 


Hel 


Hartter 


H#2 


<J/a  converter 

(•0  nr) 


1te3 


left 
t ralll no- eoge 
nap  actuator 


t«4 


etel 


rate  ltmter 


Si 


»  N/* 


rtfl 


Hotter 


rvil 


Pte2 


d/a  converter 
(80  nz) 


rte3 


traf 1 fng-eog* 
flap  actuator 


ltedef 


rtedef 
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APPENDIX  B:   MODEL  COMPUTER  PROGRAM 


*  *************************************************** 

*  **  ** 

*  **      F-18   FLIGHT  CONTROL  LAW  SIMULATION        ** 

*  **  ** 

*  **************>4****  ************  ******************** 
* 

*  *************************************************** 

*  *****INTEGRATING,IAG  ,AND     NOTCH    DIGITAL    FILTERS***** 

*  *************************************************** 
* 

MACBC    FCOT=ZINT  (FIN , Ki ,KE , IMP, FOOTZ1) 

PROCEDURAL 

IF  (IBP.NE.     1-0)     GO    TO     10 

IF  (KEEP. NE.  1.0)     GC    TO     10 

IF (TIME.  EQ.0.0)      GO    TO     10 

FOUI  =  KA*FIN    ♦     KE*FCCTZ1 

FOUTZ1=FOUT 
10       CONTINUE 
ENDBAC 
MACRO    FOUT=ZLAG  (FI  N  ,  K  A  ,  KB  ,  KC,  I  ME,  FI  KZ1  ,FOOTZ1) 

ERCCEDUR AL 

IF(IMP.NE.     1.0)     GO    TC     10 

IF  (KEEP.  NE.  1.0       GC    TO     10 

IF  (TIME.  EQ. CO)     GC    TC     10 

FOUT=KA*EIN    -    KE*FINZ1    ♦    KC*FOOTZ1 

FOUTZ1=FOUT 

FINZ1=FIN 
10       CONTINOS 
ENDMAC 
MACRO    FODT  =  Z NOTCH (FIN,KA,KE,KC,KD,KE,IMF,FINZ1,FINZ2,FOUTZ1,  FOOT  Z 2) 

ERCCEDORAL 

IFflHP.NE.     1.0)     GC    TO    10 

IF  (KEEP.  NE.  1.0)      GC    TO     10 

IF  (TIME.  EQ  .CO       GO    TC     10 

FOUT=KA*EIN    +    KE*FINZ1    +    KC*FINZ2    -    KD*FO0TZ1    -    KS*FO0TZ2 

FOUTZ2    =    FCCTZ1 

FOOTZ1    =     FCDT 

FINZ2    =    FINZ1 

FINZ1   =    Fit; 
10       CONTINUE 
ENDBAC 
* 

*  ***************************** 

*  *****FFECOENCY     2V  IB  AGERS  ***** 

*  ******  ******.*  ***************  * 
* 

MACRO       Z40=AV2CU0(Z2C,IKF) 
PROCEDURAL 

IF     (KEFP.SE.  1.0)     GO    TO    20 

IF    JTIME.NE.O.C)     GO    TO    5 

Z20Z1=Z20 

Z4CZ1=Z2C 

EEL=0. 0 

GO  TO  10 
5  IF     (IMF.  EC-  1-0)     GC    TO    10 

ZyC=ZU0Z1+DEL 

GC  TO  15 
10  ZU0=Z20 

DEL=  (Z20-Z20Z  1)/2.0 

Z2CZ1=Z2C 

ZU0Z1=ZU0 
15  CONTINUE 

20  CONTINUE 

ENDMAC 

BACBC       Z80=AV4080 (ZUO, IBF) 
PROCEDURAL 

IF     (KEEP.NE.  1  .0)      GO    TO    20 

IF     (TIME. NE. 0.0)     GC    TO    5 

ZU0Z1=Z40 

zeczi=zuc 

DEL=0. 0 

GO  TO  10 
5  IF     (IMP.  EC  1.0)     GC    TC    10 

Z£C=Z80Z1+DEL 

GC  TO  15 
10  Z8Q=Z40 
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DEL=(ZUC-ZU0Z1)/2.0 

Z40Z1=ZU0 

ZECZ1=Z8C 
15  CONTINUE 

20  CONTINUE 

E  ND  H  AC 

MACHO    Z80=AV208C (Z2C,IME) 
tfiCCEDUEAL 

IF     (KEEP.NE.  1  .C)     GC    TO    20 

IF    jTIHE.NE.  0.0)     GO    TO    5 

Z2CZ1=Z2C 

Z80Z1=Z20 

GO    TO    10 
5  IP     (IMP.  EC-  1-0)     GC    TO    10 

Z8C=Z80Z1+DEL2 

GC    TO     15 
10  Z60=Z20 

DEI2=(Z20-Z20Z  1)  /U.  0 
15  CONTINUE 

Z20Z1=Z20 

zeczi=z8C 

20  CONTINUE 

ENDMAC 

♦ 

*  ******  *******  **************** 

*  *******    RATE    LIMITEE    ******** 

*  ***************************** 
* 

BACEC    ECUT=RATIM(RIN,IME,INC) 
PROCEDURAL 

IF  (KEEP.  NE.  1)     GC    TO    100 

IF     (TIMS. NE. 0.0)     GO    TO    80 

EOUTZ1=RIN 
80       CONTINUE 

CEL1=BIN-RCUTZ 1 

DELLI»=LIMIT(-INC  .INC,  DEL 

ECUTA=ROUTZ1+DELIIM 

IF  (IMP.NE.  1.0)     GC    TC    90 

SCUTZ1=ROUTA 
90       CONTINUE 

EOUT=ZHOLD  (IMP,  ECUT2) 
100    CONTINUE 
ENDMAC 


1) 


INITIAL 
* 

*  *********************** 

*  DIGITAL    FILTER    INITIAL 

*  *********************** 

* 
INCCN  FB8Z1=. 0,?R8Z2=.0,FF9Z 
RG3Z1=.0,RG3Z2=.  0,SG4Z1=.0 
RG5Z1=.0/RG5Z2=.C.RG5F21=. 
NZAFZ1=.3  ,  NZ4Z1=. 0.NZ5Z 1  =  . 
ST16Z1=.0,ST16Z2=.O,IZ2Z1= 
NY6Z1=.0,7G6Z1=.0,YG7Z1=.0 
YV3AZ2=.0,YV3FZ1=.C,SV3Z1= 


*********************** 
PS,QC, PK,RK, YK    AND    RATS 


*********************** 

E 

CONSTANT  ?S=1761. O.CC=600.0, 
HK1=1.  177  ,RK2=3.  2*,SK7=50. 
IK7».51271.XK8=.  0  4  677. YK15 
EATE20=0. 05,EATE40=0.0  25,R 


********** 

CONDITIONS 
**'**  ****** 

1=.0,PE9Z2=.0  ,PG3Z1=. 0. ?G4Z1=.0 
,  3G4Z2=.0,?.G4FZ1=.0f3G4FZ2=.  0,.  , 
0  ,RG5FZ2=.0,RG6Z1=.0,SG6Z2=. 0... 
0,ST7Z1=.0.ST15Z1=.0,ST15Z2=.0, 
. 0,TE2Z1=.0.RSR2Z1=.0,RSR3Z1=.0 
, Y73BZ  1=.0,YV33Z2=.0, YV3AZ1=.Q, 
.0 

**  ******** 

CONSTANTS 
**  **  ****** 

?K9=-1  .  15  43,? K1 0=0. 46  47, PK 12  =  0. 
0,EK  10=1.  33,  R  K  1 1  =  1  90.  0,RK13=1  .0 
=1.01266fYK16=.0  24  69,... 
ATE80=0.0125 


5654, 
r  •  •  • 


DYNAMIC 

* 

SC£T 


**************** 

*****XNPUTS***** 
**************** 

PP1=00.0*STEP  (CO) 
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PS  1=0.  C*STEP  (CO  -00.0*  STEP  (2.  C) 
EY1=2Q .0<STEP  (C.0)-4  0.0*STEP  (2.0) 
EG1=00. 0*STEP  (0. 0) 

RG1=00.0*RAMP<C.C)-00.0*RAMP(2.0) 
YG1=0.  5*5AHP  (CO)  -  1. 5*R AMP  (2 . C) 
A  A  1  =  7.  0*STE?  (CO ) +0.0*RAHS(0.  0) 
NZ1  =  0. C*STEP  (0.0 

ny  1=0.  0*eamp  (co)  -0.0*ramp(2.  c) 

ri=qc/ps 

pic=limit  (0.  0,  1800.0, ps) 

qkf=limit  (200.c,2c00.0,qc) 

************* ************** 

*****IMPCLSE    FCNCTIONS***** 
*************************** 

ime20=impuls  (0.0,eati20) 
im?40=impols  (O.Q,5ATE40j 
i mp 80=1 m pels  (c0,eate80) 

************************************* 
*****ANGLE    OF     ATTACK    SENSOR    PATH***** 

*  4**  J************:***:********  ********** 

AA2=CMPXEL(0.  O.CO,0.74,209.0,AA1*43681.0) 

A  A 3= ZH OLD  (TMP40,  A  A  2) 

AL?HAT=0.59*A A3+1  .9 

AA6=  ALPHAT-22.0 

AA7=LIMIT  (0.  0.10CCCO.AA6) 

EK17=ZHOII(IME20,E37(H2A)*F32A(QKF)) 

P.V1=PK17*AA7 

EV1A=ZHCLE(IMP20,EV1) 

******    44****444**44*   4***4 

*****RATE    GYRC    PATHS***** 

************************* 


EG2=CM 
PG3=2H 
PK  1 1  =  1 
E2A=1. 
P2B=(1 
E2C=1. 
?G4=ZL 
EV4=LI 
PG5=PV 
EV4A=Z 
RG2=CM 
RG3=ZH 
RG4=ZN 

BG5=ZN 

BG6=ZN 

RG6  A=D 
RK6=F4 
RG6B=L 
RG7=RG 
YG2=CM 
YG3=ZH 
CSAOAT 
YG4=YG 
SNAOAT 
YG5=YG 
YG6=RA 
Y3A=(Y 
Y3E=YK 
Y3C=1. 
YG7=ZL 
YG8=RA 


PXEIjO. 
OLD  (ISP 
•65*F22 
0+EK11* 
.Q+PK1 1 
0-EK12 
AG  (PG3. 
MIT  (-80 
4*  (  (F68 

HOicria 

PXELJO. 
OLE  (IMP 
OTCH(RG 

RG3 
OTCB(RG 

RGU 
OTCB(RG 

EG  5 
EAESP(- 

IMIT  (-' 
6*RK6+5 


O.C.0,0 
8<5,EG2) 

i^O-EK 
)*(1.0- 

E2A.E2E 
.0,  120. 
(CC) *F3 
E*0,E74 
0 .0.0,0 
8C,EG2) 
3,.  165  1 
Z1,HG3Z 
U, .7589 
F21  .EGU 
5,. 6122 
F21 ,EG5 
HV16,RV 


.89,78.5,EG1*6162.25) 

EK12) 

.P2C,IMP80,PG3Z1  ,PGUZ1) 

0,?G4) 

2A  (QKF)  )  +  F40  (RIr  ES,CC,PIQ)  ) 

.8,90.0,RG1*8100.0) 

8,.  3  303  6,  .16 5 18, -.5 00 8 U, .16157, IMP 

2,RGUZ1  .RG422) 

2,-  1.  24  8  3  1,.  7  053  U,  -1.248  31  ,.4  64  26, 

F22,RG5Z1  ,EG5Z2) 

3,-.  275  0  5,-5  775  1 ,-.275  05,. 18974,  IM 

FZ2  .RG6Z1  ,EG622) 

16,  RG6) 


80,.  .  . 

IMP80  ,. 
P80, 


1 

S 

PXELJO. 
OLE  (IMP 
=  1  .C-0. 
3*CSAOA 
=(1.0-0 
4-SNAOA 
TLK  <YG5 
K15+(1  . 
15*  (1.0 
0-YK16 
AG  (YG6, 
TLM  (YG7 


CO,  10.! 

G6A+EG6! 

0.0.0,0. 

40, YG*) 

5*<  (ALPH AT* 0.0 17 4533)** 2.0) 

T 

.  12  5*  (  (ALPHAT*0. 0174533 )** 2.0) )  *ALPHAT*0. 017  4  533 

T*ZHCLG  - 

.IMP40, 

0-1  .987' 

-YK16) 


0,RG6)*LIMIT(0.0,0.12,.12-RK6) 
iE 
I. 89,  200.0,40000.  0*YG1) 


(IMP40,  RG2) 


15.  0*HATE40) 
~*YK15)  *YK16) 


Y3A.Y3E 
,IMP40, 


,Y3C,  IMPUC.YG6Z1  ,YG7Z1) 

15.0*RATE40) 


*********************************** 

*****NCRKAL    ACCELEROMETER    PATE***** 
*********************************** 
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NZ2=CMPXEL(0. 0,  C. 0. 0. 89, 2  00. 0,NZ 1*<i 0000.0) 
NZ3  =  ZriOLE  <IMP40,MZ2) 
N  Z  A = S  Z  3 

NZAR1=-50.0*LIHIT (O.Q,«.0,HZA) +280.0 
RV16=AMIN1  (NZAR  1.EK1  if 
NZ4=!IZ3-  (EG3**2)  *  (6.  8529 E -6) 
NZUA=ZHOLD(IHP20,BZU) 

HZAF=ZINT  (NZ 4  A,  0.25,  0.75,  IMP2 0  ,NZ AFZ  1) 
P5A=  (1.  0+EK9*  (  1.0-PK  10)  ) 
Pf3=(1  .0  +  PK9)  *  (1  .C-FK10 
F5C=1.0-EK10 

NZ5=ZLAG  (NZ4.E5A  .E5E,P5C ,IMP4 0  ,NZUZ1  rNZ5Z1) 
EV2=LIMIT  (-1(5.  0,  10.0, :iZ 5) 
NZ7=PV2*3.5*F32A  (CKF) 
EV2A=ZHOLE(i:iE20,EV2} 
* 

*  *************  *******  ******* ******* 

*  ****LATEEAL    ACCELEEGMETSR    PATE**** 

*  ********************************** 

* 


ZHOLD  (IHPUO,  RG3) 


* 

*  ******************************** 

*  ********£ ilot    INEOT    P&IH ******** 

*  ******************************** 
* 

PP2=DEADSF(-2  .C.2.0.EP1) 

EE3=C.1FXEL    0.  C,C.0.0.14  .27.3  , EE2*106  .47) 

EE4=RS ALEL(0. 0  .  7. S36 5E- 3 , PP 3) 

PP5  =  ZHOLE  (IflP8C,EE4) 

EE6=P?5*  (7.0*  (0.2*AES(?P5))) 

PP7=LIMIT  (-25.C,5C.C,PP6) 

EV3=P?7 

EV3A=ZHOLC  (IME20.EV3) 

PR2=DEAD£E  (- 2  .  C ,2  .  0  .  ER1 ) 

PE3=C.1PXEL(0.  0,  C.O,  0.7,  24 .3 , PE2* 1 61 .  OH) 

PR4=LIMI1  (-3.0,3.0.  FR3) 

EE 5= REAL  EL (0.0,0. 0079, PR4) 


EE5=REALEL(0.U,0.G07  9,PR4) 

PR6=ZH0LE  (IMP8C.EE5) 

PR7=DEADSF  (-0.C25  ,0.  C25.PE6) 

EE8=PR7*  (AES  (EE7)  *EK1+SK2) 

?R9=ZNOTCH(PR8  , . 6555 .-1 . 5 51 2  ,  .  8234 , -  1 . 551 2 , . 6  789, IMP80 ,. 

PR8Z1,EE8Z2\PR9Z1,PBSZ2) 
PY2=DEADSP(-7.0,7.0,EY1) 

FY3=CMPXEL(0.  0,0.0,  1.5,  66.6,4  4  25.6*PY2) 
?Y4=ZH0LE  (I3PUC, E13) 
EY5=DEADSE (-1 .  0,1 .0, PY4) 
YK9=0. 00072*F 10 (QC) 
YK10=0.234*F1  C  (QC) 
EY6  =  PY5*  (YK9*AE£  (EY5)  +YK10) 
YV1=PY6 

EY7=PY6*  (0.5-  (E17 fALFHAT)  *F1  U(HI))  ) 
PY8=LIMI1  (-3  0.C,3C.0,PY7 


* 

*  ************************* 

*  *****STAEILATCE    PATH***** 

*  ************************* 

* 

EK19=ZHOIE(IME20,E12  (RI,  PS)  *  F32A  (QKF)  *3  .5) 

PK21  =  ZHOLB(i;iE20,r  12  (RI,PS|  *F;2A  (QKF)  *F68  (QC)  ) 

PK22=ZHOLE(iaE20,F12  (RI  ,PS)*F32A  (QKF)  ) 

ST1=PV3A*EK22 

S12=(PVUA*PK21)-S11 

ST  3=  (PK19*PV2AJ  +  PV1A+ST2 

S1U=LIMIT  (-10CCC.C.0.0,ST3) 

ST5=LIMIT  (0.  0,  10  0  00.  0,5  T  3) 

ST6=(ST4*F23  (  AIPH AT)  ) +ST5 

ST  8  =  1 3  PL  (0.0.  0.0  5,  ST  7) 

S16A=ST€-C.5*ST8 

ST7=ZIMT  (ST 6  A, 0.05, 1.0, IHP20, ST7Z1) 
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ST9=A7 
ST10=N 
SI11=S 
ST12A= 
ST12=Z 
ST  13= A 
ST14=S 
S115=L 
ST16=2 
I 
ST17=S 
ST18=L 
DS1=RV 
DS2=LI 
D£3=DS 
ST2R1= 
SI2R2= 
CS4=LI 
DS5=RA 
R£T1=S 
RST2=L 
R£13=Q 
RSTDEF 
LSI1=S 
LST2=L 
L£I3=Q 
LSTDEF 


20  4C  (LIMIT  (-5  0.0,2  5  .0,S*I8)  ,  IMP  20) 

Z7+ST9 

11 C+PV1 

ST11+((EG:-*YG3l*F107(QC)) 

HOLD  (IMF40,ST12A) 

V4060  (ST12, IMP40 

T1  3  +  DG5 

IMIT  (-25.0,25  .0.ST1  H) 

NOTCH (ST  15,. 69 0  84, -.9 90 6 €,.66 3 12, -.99 068,. 353 96, 

MP8C,ST15Z1  .ST15Z2.ST16Z1.ST16  12) 

T16-(P73*F3*A  (CKF)  ) 

IMIT(-24.0,10.5,ST17) 

2+EV1 A*1 .33 

MIT  (-50.0,50.0,DS1) 

2*R*7 

17.25-AES (ST2+6.75) 

LIBIT(-5C.0,F41  (A  LP  HAT,  EI)  *10.  0,ST2R1) 

MIT  <-ST2F2,ST2E2,  DS3) 

TLE(DSU,IMF60,35.0*EATE8C) 

T1 8-DS5 

IMIT(-24.0.  10.5.RST1) 

NTZR  (RATE80,RST2) 

=CKEXFL (0.0,0.0,0.7,40. 0,RST3*  16  0  0.0) 

11 8+DS5 

IMIT  (-24.0,  10.  5.LST  1) 


NTZR  (RATE36,LSt£) 
=CKFXPL  (0.0,0. 0,0, 


7,40. 0,LST3* 1600.0) 


****AILE5CN    PATH**** 
***** ********  ******* 


EV2A=P 
RV2=LI 
RV1  1=F 
EK17=F 
RV7=AM 
Y  V  1  F.=  L 
RV1=LI 
RV1A  =  A 
AL1=RV 
AL2=LI 
AL3=AL 
AL4=LI 
AL5=DE 
EAI1=L 
RAL2=R 
RAL3=Q 
RALDEF 
LAL1=L 
LAL2=R 
LAL3=Q 
LALDEF 


R9*  ((F 

MIT  (-5 
31  (F.I, 
6(FI,? 
IN  1  (RK 
IMIT  (0 
MIT  (-2 
V4C5Q ( 
2+E71A 


2*F35 


MIT  (-5 

MIT  (-5 
ADS?  (- 
IMIT(- 
AIIM(R 
NTZE  (R 
=CMFXP 
IMII(- 
ATLM(L 
NTZR  (R 
=CSEXP 


p  e\  *F 

S  7  E  V  1 
1  7*F3 
.0,1. 
00.0. 
HV1.I 
*  1 .  33 
0.C.5 
ALPHA 
O.C*F 
42.0, 
42.0, 
AL1.I 
ATE80 
L  (0.0 
25.0, 
AL1.I 
AIE80 
L  (0.0 


ES)  +F13  (RI,  P£))*F7  (CC,PS)  )-RG7 

C.0,RV2A) 

34  (ALPHAT) 

5  (ALPHAT)  ,  RK17-F101   (RI  ,  PS  ,  NZAF)  ) 
0,F29(AL?HAT) )*YV1 
200.0,  YV1R) 
2E40) 


C.0.AL1 ) 

T)  *E36  (RI,QC,FS) 

41  (ALPHAT, SI)  ,50.0*F4  1  (AL 
42.0,AL4)+AL4 
25.0,AL5) 

2P8C.1  00.0*RATE80) 

,EAL2) 

,O.C,0.7.40.C,-160  0.0*EAL3) 

42  .  C   AL5  ) 

«P80|100.0*RATE30) 
,LAL2) 

,0. C,0. 7,40.0, 16  00. 0*LAL3) 


PHAT,RI) , AL3) *0. 5 


******************* 

****RUDDER    PATH**** 
******************* 


ESE1 
RSR2 
Y5A  = 
Y5E= 
Y5C  = 
RSS3 
SSR4 
YV2= 
YV3D 
IV3C 
YV3B 
YV3A 

Y73  = 
EC1= 

RC2  = 
EE3= 
RRD1 


=  LAL 
=  ZHO 
(1.0 
(YK7 
1.0- 
=  ZLA 
=RSR 
F4  2  ( 
=ZHO 
=  YG8 
=  LIM 
=ZNO 


2  +  RAL 
LD(IM 
+  YK7* 
♦  1.0) 
YK8 
G  (RSR 
3*F30 
BSR4) 
LD(IM 
*F96  ( 
IT  (-3 

(Y73A 

YV3-P 


2  +  C 
E40 
(1. 
*(1 

2,Y 
(CC 


Sf*2 
.ESP. 
0-YK 
.Q-Y 

52, Y 
,PS) 


tc.ai: 


ZLAG 
Y72  + 
RD  1* 

LIMI 
=LIMIT(-3 


51^8 


P40 

C6 

V3E 
3EZ 

/*" 
Y8 

C(5P 
c.5 


EG  3 


,.13 
1  ,YV 

4837 

30.0 
,30. 


53, Y5C, IMP40,RSR  2Z1 ,RSR3Z1) 

*F38  (ALPHAT, RI, PS) 

)  *ZHCLD(I£E4  0,?G3)*F108(QC) 

NH- LIMIT  (-30.0,  30.0,YV3D) 

0,Y73Q 

876 ,.277  5  2..1387  6t-.7  914,.346  42,IMP40, 

2BZ2,  YV3AZ1.Y73AZ2) 

,-.44  33  7,.10326,IMF40,YV3FZ1,YV3Z1) 


,RD2) 
C,RD3) 
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RRD2=A  V4CS0  (RRD1  .IMPUO) 
HBE3=QNTZE(H4TI80,BBE2) 

RRDDEF=CHPXPL fC. 0,0. 0,0.7, 4 O.0, 16 00. 0*RBD3) 
LED1=LIMIT (-3  0.0,30. C.RD3) 
LHD2=AV4CE0(LED1,IHE40J 
LEE3=QNTZE(RATE8Q,LRE2[ 

LEEDEF=CKPXPL (0.0,0.0,0. 7,40. C, 16  00. 0*LRD3) 
* 

*  *****LEAEING    EEGE    FLAP    PATH***** 

*  ******************************** 

* 

LE2=ZINT  (ALPHAT.0.06  25,0. 937  5 ,IHP40 , IE2Z1) 

LE2T=ZHOLE (IMP20,LE2) 

LE2T1=F27  (LZ2T.RI) 

LE2T2=LIMIT(O.0.F28  (CC)  . LE2T1 ) 

LE3=LIMIT  (0.  0,E29U(RI)  ,  LE2T2) 

LS4=RATLM  (LE3.IHP20. 18.0*RATE20) 

LE5=AV2080  (LE4.IME20) 

DLE1=RV2*F93 (HI ,PS , NZAF) 

DLZ2P1=LE5+3.0 

DLE2=DEAES?(-EIE2r1,  ELE2P1.DLE1) 

DLZ3=DLE1+DLE2 

RLE1=LZ5+ELS3 

R LE 2=LIMI1 (-3.0,3 3.0 -RLE  1 ) 

RLE3=RATIH(RLZ2,IH?eo.18.0*RAlE80) 

RLE4=LIMIT    -3.C.33.0  ,RLE3) 

ELE5=QNTZE  (R  AT  E80  ,  RL  14) 

RLZDEF=CMPX?L  (0. 0, 0. 0,0. 7, U0.C,RLZ5*  1600.0) 

LLE1=LE5-DLZ3 

LLE2=LIMIT(-3.0,33.0,LLS1) 

LLE3=R A TIM (LIE  2 ,1  MP  8 C , 1 8 . 0* RATE8  0) 

LLE4=LIMIT(-3. 0, 2  3.0.LLZ3) 

LLZ5=QNTZF  (R  ATE80  ,LIZ4) 

LLEDEF=CMPX?L  (0.0,0. C ,0  .  7 ,40 . 0  ,LLZ5*  1600.0) 

*  ********************************* 

*  *****TRAILING    EDGE    FLAP     PATH***** 

*  ********************************* 

* 

TE2=ZINT  (ALP HAT- 0.0 3 125, 0.96  3  75, IMP4  0.TE2Z1) 
TE2T=ZHOLD(IHP20,TE2) 

TE2T1=F24 (RI.1E2T) 

TE3=LIMIT  (0.0,E25(CC),TE2T1) 

TE4=AV20eC(TE3  .IKE20) 

DTE1=RATLH(RV2*RV11,IMP80, 18.  0*RATE8  0) 

R1E1=TE4-DTE1 

ETE2=LIMIT  (-8.  0,45.0,  RTZ1) 

RIE3=QNTZR  (R  ATZ3C  ,RTE2) 

BTEDEF=C2PX?L  (0.0,0.0,3. 7,40. C,RTZ3*  16C0.O) 

LIE1=TS4+DTE1 

LTE2=LIMIT(-8.  0r45.0,LTE1) 

LIE3=QNTZR  (R  ATE80  ,LIE2) 

LTEDEF=CMEXPL  (0.0,0.0,3.7,40. C,LTE3*  160  0.0) 
TERHIN  AL 
METHOD    EKSFX 

TIMES    FINTIM=4.00,?RDEL=0.25.OOTDEL=0.25,DELT=0.0125 
TITLE    RUEEER    RESPONSE    TO     +/-20LE    YAH 
PRINT    RRDDEF,LEDDZF 
PAGE    XY  PLCT 

*  OUTPUT    TIKE.  RLEDEF  (-5.0,5.  0),  LIEDEF  (- 5  .  0,5.0) 

*  LAEEL    OIFF.L.  Z.     FLAP     RESPONSE    TO    +/- 1  OLB    ROLL,     'WITH    RG    FEEDBACK 

*  LAEEL  »=1.40,H=500CFT 

*  OUTPUT  TIME,  RTZDZF  (-5.0,5.  0)  ,  LTEDEF  (- 5.  0,5.0) 

*  LAEEL  DIFF.T.E.  FLAP  RESPONSE  TC  +/- 1 OLB  ROLL,  WITH  RG  FEEDBACK 

*  LABEL    M=1  .4O,H  =  500C    FT 

OUTPUT    TIME.REEEEF  (-5.0,  5.  0)  .  LBEDEF(-  5.0,5.0) 

LAEEL    RUDDER    RESPONSE    TO    +/-20LE    YAW,     WITH    YG,NY,RG    FEEDBACK 
LAEEL    M=0.81, H=5000ET 
END 
STCP 

FUNCTION    F4(QC,PS) 

EEAL    LIMIT 

F4T1=0.36-0.  0C12*CC 

F4T2=LIMIT(2O0.,1C0C.,PS) 

F4T3=0.2-C.00  0  3*F4T2 
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F4T4=0.  1-0.000  1*F4T2 

IF(.NOI.F4T2.LE.5C0.  )    GO    TO    10 
F4T5=F4T3 
GO    TO    2C 
10  CONTINUE 

F4T5=F4T4 
20  CONTINUE 

F4=LIMIT  (FUT5,C.3,?UT1) 

5ETURN 

END 

FUNCTION    F6  (RI,PS,RV11) 

REAL    LIMIT 

F6T1=0 .2  2*LIMIT  (0.0,  1.  12,  31)  +  ( 1.0 HE-  4)  *FS 

F6T2=LIM IT (0.0. 0.43  1  .F6T  1 )  -0  .  25*RV1  1 

F6=.5  1 1 -LIMIT  (C. 311,0. 431,  F6T2) 

RETURN 

END 

FUNCTION    F7(QC,PS) 

REAL    LIMIT 

F7T1=0.0  5*LIMIT(25.0, 125 . 0, QC) +3. 75 

F7T2=QC-325.0 

F7T3=PS-628.0 

F7T4= (LIMIT (-3  25. ,3  2  5.,F7T2)  *IIMIT (- 386 . , 0. , F7T3)  ) *  (2. 7 1E- 5)  +10.0 

IF(. NOT.  CC.LE. 325.0)     GO    TO    10 
?7=F7T1 
GO    TO    15 
10  CONTINUE 

F7=F7T4 
15  CONTINUE 

RETURN 

END 

FUNCTION    F10(QC) 

REAL    LIMIT 

F1CT1=1 .C-O.00C425*CC 

F 10=  LI  MIT  (0.  35,1.  Cr  7 1011) 

RETURN 

ENE 

FUNCTION    F12  (SI, PS) 

HEAL    LIMIT 

F 1211=81**2*9.625-. 025* RI+1 .0 

F12T2=PS*7. 96  9E-4+0. S4 

F12MAX=LIMIT  (1  .  C  ,  8.  0  ,F1  2  T2) 

F12T3=LIMIT  (1.  0,  F  12  MAX,  F1  2T  1  ) 

F12T5=LIMIT  (.5,1.  35.  EI) 

F12T6=F12T5* (0 . 00S52 *PS+ 4  .0  4)     +  (-PS* 0. 0  039  6- 1 .  18) 

F12T4= LIMIT (1.C,8.0,F12T6) 

IF  (.NOT.  SI.  GT.  0.5)    GOTO    5 

F  12    =    F12T4 

GO    TO    6 

5  CONTINUE 
F12    =    F12T3 

6  CONTINUE 
RETURN 
ENE 

FUNCTION    F13  (HI, PS) 

HEAL    LIMIT 

?S13L=LIMIT(2CC.C,21 16.0, PS) 

F13T1=  (445.  0* LIMIT  (242.  0,  628.  0,PS13L)  )  *  (2.  39 E- 4) 

IF  (. NOT. PS13L. IE. 800.0)      GO    TO    10 
F13T2A=0.054+psi3L*(3.59S-4) 
F13T2E=~0.  T/45-PS  13L*  (  S.4E-5) 
GO    TO    20 
10  CONTINUE 

F12T2A=O.C0122*PS13L-0.6  37 

F13T23  =  0.  192-PS13L*  (5.4  2E-4) 
20  CONTINUE 

F13T2C=0.  152  +  PS13L*(5.34E-5) 

RI13L=LIMIT(0.C.2.8,EI) 

F13T3=  F13T2C+SI13L*  (  F13 T2 E+RI 13L*F1  3 T2A) 

F  13=  LI  MIT  (0.  13,F12T1,F13T3) 

RETURN 

ENE 

FUNCTION    F17  (AIPHAT) 

EEAL    LIMIT 

F17T1=ALPHAT-13.0 

F17=LIMIT  (0.  0,4.0, F17T1) 
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BETURN 

END 

FUNCTION  F22(CC) 

REAL  LIMIT 

F22=0.0  167*  (-9  00.C+LIMIT  (600.  C,900.0,QC)) 

RETURN 

ENE 

FUNCTION  F23  (ALPHAT) 

SEAL    LIMIT 

F23  =  3.  14  2  5-0. 142  9*LISIT (15.0,21.998, ALE  HAT) 

RETURN 

END 

FUNCTION    F24  (HI,  ALPHAT) 

REAL    LIMIT 

F24L1=  22.538-20.5  1*LIMIT(0.  27,  0.66,  R  I) 

F24L2=32.76-3  6.0*IIMIT(0. 66 ,0  .91 ,RI) 

F24T1=LIMIT  (0.  C,1C00  0.0,  ALPHAT) 

F24T2=  ALPHAT-  (  14 .  67  6  9-7  .  6  92  3*LIMIT  (  0  .  27  ,  0 .  9  1  ,RI)  ) 

F2  4T3=-2.C*LIMIT  (C.  C  ,10  00  0.  0  ,  F2UT2) 

F24T4=1.4*(F24T1  +  E2  4T3) 

IF    f. NOT. HI. GT. 0.66)    GO    TO    10 
F24L=F24L2 

GO    TO     15 
10  CONTINUE 

F24L=F24L1 
15  CONTINUE 

F24=LIMIT  (0.0,F2  4I,F24T4) 

RETURN 

END 

FUNCTION    F25(CC) 

REAL    LIMIT 

F 25=4 7 .6  26-0.0  5106*LIMIT ( 600. 0 ,835. 0 ,QC) 

RETURN 

END 

FUNCTION    F27  (ALPHAT, EI) 

REAL    LIMIT 

F27=1. 3  26*(ALPHAT+7.8584-17 .8  6*LIMIT  (0. 44 , 0.  63 , RI)  ) 

RETURN 

END 

FUNCTION    F28  (QC) 

REAL    LIMIT 

F28=44.5  5  1-0.0  40  5  6*LIMIT ( 26  0. C, 95  0. 0 ,QC) 

RETURN 

END 

FUNCTION    F29U  (RI) 

REAL    LIMIT 

F  29  0=8  7.  2  82  5- 7  6.  25*  LIMIT  (0.7 , 1 .146, RI) 

BETURN 

END 

FUNCTION    F30  (C.CPS) 

REAL    LIMIT 

F30T1= 1. 0  803-  (LIMIT (  140. 0,8 00.  G,QC) )  *0. 00 12878 

F3  0T2=1  .  7C35-  (LIMIT  (242.  0,4  9  8.0,  PS)      *0.  0  0  33203 

F30  =  LIMIT  <F3  0T2,0.9,F30T1  j 

RETURN 

END 

FUNCTION    F31  (EI, PS) 

REAL    LIMIT 

F31T1=  (LIMIT  (0.0,  2  116.0,  PS)  )  *  (6.536E-5) 

F  2  1T2=  1.449-]  LIMIT  (1455.  0  ,21  16.0,  PS)  )  *(5.  1  E-4) 

IF  (.NOT.  BI.LE.  C.  717)     GO    TO    10 
F31T3=F31T2 
GO    TO    20 
10  CONTINUE 

F31T3=0. 112 
20  CONTINUE 

F21T4=RI-C717 

F31T5=  (LIMIT  (-0.7  17,0.662,?31T4))*F3  1T3 

F31T6=0.2  24  7-F31T1-?31T5 

F31=LIMIT  (0.  0,0.  16,F21T6) 

REIURN 

END 

FUNCTION    F32A  (CKF) 

F32A=100.0/QKF 

RETURN 

END 
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FUN'CTICN    F34  (  ALPHAT) 

REAL    LIMIT 

F3UT1=  1.5-0.  2*  AES  (ALFHAT-2.  5) 

F3  4=LIMII  (0.0,1.0,F3  4T1) 

EETURN 

END 

FUNCTION    F35(ALPHAT) 

REAL    LI  KIT 

F35T1=0. 955-0. C3 25* ALPH AT 

F35T2=C.5*0.0  4  221*ALEHAT 

IF  (.MOT.  A1PHAT.LE.6.  C)     GO    TO     1C 
F35T3=F35T2 
GO    TO    20 
10  CONTINUE 

F  "'CT3  =  1?3  CrT1 
20  CONTINUE"' 

F35  =  2.C*LIMIT(C.175,0.5,F35T3) 

RETURN 

END 

FUNCTION    F36  (RI,QC,FS) 

REAL    LIMIT 


F36T1=  LIMIT  fO.  0,0.8,  RI)  +  PS*  (  1 .  69E-4)  -.7  68 
F36T2=  (4  4C.0-LIMIT  (4  40.  0,640.0,PS))*7.46E-3 
F36T3=  (LIMIT (0.7, 1.9,  SI)  )  *54 0 . C  +  955 . 0-QC 


'  (44C.0-LIMIT  (440.  0, 

■  (LIMIT (0.7, 1.9,  SI)  ) 

JLIMITjO.O, 200.0, F3 


,9,31) )*5  4  0. C  +  955. 0-QC 
F36T4= JLIMIT [0 . 0  ,  20  0 . 0, F3  6T3)  ) +0.005 
F36T5=F36T1*F36Tz+1 .C 
F36T6= LIMIT (0.0,1.0,F36T5) 
F36=AMIN1  (F36T4,?36T6) 
RETURN 
END 

FUNCTION    F37(NZ2) 
REAL    LIMIT 

F37  =  2.  5-0.5*LIMIT  (3 . C,5 . 0 ,NZA) 
RETURN 
END 

FUNCTION    F38  ( ALPH AT , RI, PS ) 
REAL    LIMIT 

FSSI =L I MIT(0. 225, C. 549,31) 
IF  (.NOT.  ALPHAT.LT.  16.0)      GO    TO     10 

FS=-0.243C8*FSEI+0.03  9  67 

AOAM=-10.  0 

AOAB=-0.45 

GO    TO    20 
10  CONTINUE 

FS  =  -0.  13846*FSEI-0.0 1774 
ACAM=- 5.467 
ACAE=-0. 025 
20  CONTINUE 

F38T1=LIMIT (-5. 0, 20. a, AL? HAT) 
F3£T2=LIMIT(-5.C,25.0,ALPHATl 
IF  (.NOT.  EI.LE.  0.  334)     GO    TO    30 

F33T3=F38T2 

LL=-0.98 

GO    TO    40 
30  CONTINUE 

F36T3=F3ET1 
LI=-1. 2 
40  CONTINUE 

F38T4=  (F38T3+ACAM) *FS  +  AOAE 

F38T5=LIMIT(LI.0.3.F3  8T4) 

F38T6=  0. 048-0.  C4  9£*LIMIT  ( -5 . 0  ,25 . 0 , ALPH AT) 

F3  8T7=LIMIT  ( -  5  .  0  .  25  .  0 ,  ALP  HAT)  +0  .  0  1*L  IMIT  (  240 .  0  ,  3  90 .  0  ,?  S)  -3 

F3  8T3=LIMITJ0.C,2.5rF38T7)-2.5 

F38T9=RI-0.72 

F38T1 0=0. CI 6 8*11  HIT  (CO. 2  .0  8,F38T9) 

F3  8T11=F38T10*E38T8+F38T6 

IF  (.NOT.  AIPHAT.IE.10.0)     GC    TO    50 

F38=F3ST1 1 

GO    TO    60 
50  CONTINUE 

F38=F38T5 
60  CONTINUE 

RETURN 
END 

FUNCTION    F39  (  ALPHAT) 
REAL    LIMIT 
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F3  9= (LIMIT  (13 . C,25.C,ALPHAT)  )  *0. 0833 3- 1 . 08  32  9 

EETURN 

END 

FUNCTION    F40 (EI,FS,QC,PIQ) 

REAL    LIMIT 

F40T1= 3. 25-3.0 

F40T2=C. 65625- 

F40T3  =  -0. 26177 

F40T4=LIMIT  (0 


F40T5=(F4CT4*  ( 

F40T6=F40T3+  (F 

F4CT7=F4CT4* (5 

F4CT8=  (0.  16923 

IF     (.NGT.BI.GE 

F40T9=F40T8 

GO    TO    30 

20  CONTINUE 

F40T9=F40T8 
30  CONTINUE 

IF     (.  NOT.  [QC.  G 
F40T1C=F40T 
GO    TO    50 
40  CONTINUE 

F40T  10=F40T 
50  CONTINUE 

IF    7. NOT. EI 

FUQ=FUCT 

GO    TO    80 
60  CONTINUE 

IF     (.NCT.EIQ.L 
F4Q=F4CT5 
GO    TO    80 
70  CONTINUE 

F40=F40T10 
80  CONTINUE 

EETURN 

END 

FUNCTION    F41  (A 

REAL    LIMIT 

F41T1= 

F41T2=(. 23533* 

F41T4=0.  5875-  ( 

IF  (.  NOT.  AL2HAT 
F41 13*1. Q-F 
GO    TO    20 
10  CONTINUE 

F41T3=  1.  0-F41T 
20  CONTINUE 

F41  =  LIMIT  (F4  1T 

RETURN 

ENC 

FUNCTION    F42(I 

REAL    LIMIT,     IN 

FU2=(A3S  (LIMIT 
*  *LIMIT<-18. 

RETURN 

f  n  r 

FUNCTION    F45 (Q 

EEAL    LIMIT 

F45T1 

F    " 

F45  =  1 . C-FU5T1 * 

EETURN 

END 

FUNCTION    F68(Q 

REAL  LIMIT 

F66=-0.0C2977* 

EETURN 

END 

FUNCTION  F90 (E 

REAL  LIMIT 

F90T1=628.0-?S 

F90T2=RI-C.72 

F90=  (LIMIT  (0.  0 

RETURN 

ENC 


♦LIMIT 
-0C131 
+  (9.63 
C,  1C5. 

-EIS*( 

40T2*L 
•6746E 
-3.86U 
.0.75) 
*F4CT1 


(0.75 

25*LI 
5E-4) 
0,QC- 
1.174 
IMIT  ( 

-4+  1 . 

153-5 
GO    T 


,0. 

MIT 
*LI 
335 
28  E 
0.7 
98  4 
*pj 
0"  2 


85, RI) 

(500.0.  1800. O.PS) 

MIT  (500-  0,1800.0, PS) 

-6U  (1.5  238E-3)  )  )  +0.  475- (6  .5E- 4)  *PIQ 
5.0.35,3  1)) 

1E-7*FIQ)  +  (0. 169  23-3. 84 615 E-5*PIQ) 
O)+(F40T4*(1.67247E-3-9.29  152E-7*PIQ)) 


E. 440. 0. AND. RI.GE. 0.75) )     GO    TO    40 
6 


GT. 980.0)     GO    TO    6C 


E.5C0.  C)     GO    TO    70 


LPHJT,EI) 

4  1T1=  (LIMIT (8.0, 22.0. ALP  HAT)  -8.  0)*0. 052571 

4  1T2=  (.2353 3* LI  MIT (.27,. 48. RI) -. 009 97)  * (LIMIT (8 . , 22. , ALPH AT ) -8 . ) 

41T4=<3.5875-  (LIMIT  (  0.2  7  ,  0.  3  9  ,  FI)  )*1  .25 

.IS. 8.  0)     GO    TO     10 

41T1 


4, 1.0, F41T3) 


8) 


'45T1=  (LIMIT (5 
?45T2=  (LIMIT  (6 
'45=1 . C-F45T1 * 


(- 18.  3  3. 18.  33, 1  N) )  *0.  0  479 +0.76  1) 
zl,  18. 33, IN) 


C,ES) 

00.  C,1 000.0, QC) -500.0)  *0.001 

25.0, 1450. 0,?S    )*0.  00  1 212  1-0.7  57  5  8 

F45T2 


C) 
(-480.0+LIMIT  (2  60.0,4  8  0.0,CC)) 

I, PS) 

,386.0,F9  0T1)  *LIMIT(0.  C , 2. 08 ,F90T2) ) *0. 02055 1+ 1 5 . 5 
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10 

20 


10 
20 


157.0,F9  3T2) 


ENDJCB 


FUNCTION    F93 (EI,ES,NZAF) 

REAL    LIMIT, NZAF 

F93T1=  (RI+0. 0  00  18*PS-0. 69) *0. 12857 

F93T2=PS-€28.0 

F9  3T3= 1.0-0. 1*LIMIT(0.0,5.0,NZAF) 

F9  3T4=LIMIT (O.C.0.06  ,F9  3T1) 

F9  3T5=  (6.36  94E-3)  *LIMIT  (0.0, 

F9  3=F9  3T4*'?93T5*FS3T3 

BETURN 

END 

FUNCTION    F96  (QC,PS) 

REAL    LIMIT 

F96T1= 1.51-0. 0C44+LIMIT (0.0,168.0,QC) 

F96T2=CC*0.0  00  88+0. 3  6 

IF (.NOT. C.C.LE  .500.0)     GO    TC    10 

F96  13=F96T1 

GO    TO    20 
CONTINUE 
F96T3=F96T2 
CONTINUE 

F96T4= 0.000777  2*  (6  28.0-LI  MIT  {  212.  0,  6  28.  0,  ?  S)  ) 
F96T5=F96I3+F96I4 
F96=LIMIT  (0.0,  1.  1,F9  6T5) 
RETURN 
END 

FUNCTION    F101  (EI, PS, NZAF) 
REAL    LIMIT, NZAF 
F101T1 =51-0.58 
F101T2  =  LIKIT(0.C,5.0,F10  1T1) 
F1C1T3=PS-498.0 

F101T4=F101T2*LIMIT(0.0,8  35.0,F101T3) 
F101T5=7  9E.54-C.35S4  6*PS 
F101T6=LIKIT  (165.2,27  5.  5,F101T5) 
F1C1T7=(5.4UU6E-5)*LIMIT  (0.  3,F10  1T6,  F10  1T4) 
F1C1T8=AES  (NZAF) 
F101T9  =  LIMIT  (0.0  ,5. 0,F10 1T8) 
F1C1=F 101T7*F 1 C1T9 
RETURN 
END 

FUNCTION    F107  (CC) 
REAL    LIMIT 

F107T1=  (  (-5.  7  HE- 7)  *ABS(QC-750.0))  +  (8.4E-4) 
F1C7=LItfIT(O.C,10C0O.0,Fl07T1) 
RETURN 
END 

FUNCTION    F108 (CC) 
REAL    LIMIT 

F108T1=LIMIT  (5  0.  0,39  0.0,QC) 

F1Q8=(F1C£T1**2.0)*7.48E-3+0. 015  57-F  108T1*6. 43  5E-5 
RETURN 
END 

FUNCTION    F112  (EI) 
REAL    LIMIT 
IF  (RI.  LT.O.  28  -CR.  EI.  GT.  0.76)     GO    TO    10 

F112=4.0 

GO  TO  20 
CONTINUE 
F  112=0.0 
CONTIN  UE 
RETURN 
END 

FUNCTION  F113  (JLPEAT) 
REAL  LIMIT 

F  113=0.  16  667*LIMIT(12.0, 
RETURN 
END 

FUNCTION    F114  (RI) 
REAL    LIMIT 

F  111=0. 89286* LIMIT (0.16 
RETURN 
END 


18.0,  ALPHAT)  -2.  0 


0.30 ,EI) -0. 1 4286 
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